#ifndef _const_nucl_
#define _const_nucl_
const double Pi = 3.141592653589793;
// average neutron and proton masses.   2006 value of average (m_p+m_n)/2
//(/938.918725_dp, 938.918725_dp, 938.918725_dp, 1232.0_dp, 1232.0_dp/)
const double P_massave = 938.918725;

const double hbarc = 197.3269678792965; // Mev.fm
const double Mpic2 = 138.04;            // MeV  139.57;     // Mpi c2

// New LECs from Luigi */
/*
>> c1=-0.81
>> c2=3.28
>> c3=-3.0
>> c4=3.40
>> cD=1.85    cD and cE have been fitted on triton binding energy and beta
>> decay
>> cE=0.313
*//*
const double LEC_C1 = -0.81; //!! GeV^-1
const double LEC_C3 = -3.0; //!! GeV^-1
const double LEC_C4 =  3.4; //!! GeV^-1
const double LEC_CD =  1.85;
const double LEC_CE = 0.313;
//

const double LEC_CD = -1.0;
const double LEC_CE = -0.34;
const double LEC_C1 = -0.81; //!! GeV^-1
const double LEC_C3 = -3.2;  //!! GeV^-1
const double LEC_C4 = 5.4;   //!! GeV^-1
// PRL 99, 042501 (2007) When CD = - 1.0 CE = - 0.34 From figure FIG. 1


const double LEC_CD =  0.816806;
const double LEC_CE = -0.039575;
const double LEC_C1 = -1.121521; //!! GeV^-1
const double LEC_C3 = -3.925006; //!! GeV^-1
const double LEC_C4 =  3.765687; //!! GeV^-1
*/
//  NNlosat



const double Cut_LambdaX = 700;
const double Mass_pi = 0.704110550591284; // 138.94/hbarc;
const double F_pi = 92.4;

// is setted in Wave_Fun
// const int k_mesh=16;
// const int p_mesh=16;

// int nlamx,lmax,nmax;
// double hbar_omega;
// const double Cut_Lambda = 500;
class LEC_CONST {
public:
  double LEC_gA;
  double LEC_C1; // = -0.81; //!! GeV^-1
  double LEC_C3; // = -3.2; //!! GeV^-1
  double LEC_C4; // =  5.4; //!! GeV^-1
  double LEC_CD; // =  - 1.0;
  double LEC_CE; // = -0.34;
  LEC_CONST();
  LEC_CONST(double c1, double c3, double c4, double cD, double cE)
      : LEC_gA(1.29), LEC_C1(c1), LEC_C3(c3), LEC_C4(c4), LEC_CD(cD), LEC_CE(cE) {}
};
#endif
